A Superconducting Gap in an Insulator 
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We present tunneling spectroscopy and transport measurements on disordered indium oxide films 
that reveal the existence of a superconducting gap in an insulating state. Two films on both sides of 
the disorder induced superconductor to insulator transition (SIT) show the same energy gap scale 
at low temperatures. This energy gap persists up to relatively high magnetic fields and is observed 
across the magnetoresistance peak typical of disordered superconductors. The results provide useful 
information for understanding the nature of the insulating state in the disorder induced SIT. 
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Increasing the disorder of a metallic system causes the 
localization of its electronic wave function. If the ground 
state of the system is a superconductor, increasing the 
disorder leads to a transition from a superconductor to 
an insulator (SIT). This transition has gained increasing 
attention lately due to the experimental observations of 
a number of dramatic features near the SIT such as sim- 
ple activated temperature dependence of the resistance 
on the insulating side [l[, a large peak in the magneto- 
resistance peculiar I-V characteristics 0, Q and 
traces of superconductivity at temperatures above Tc 0" 
[l^ . Another reason for renewed interest in this field is 
that the SIT may be a basic realization of a quantum 
phase transition which occurs at T = as a function 
of disorder or magnetic field and is driven by quantum 
rather than thermal fiuctuations. Despite this growing 
interest, the mechanism of the SIT and, in particular, 
the nature of the insulating state are not understood. 
Recently a few indications for the presence of electronic 
pairs in the insulator have been reported [13, [3, [3], 
inspiring further theoretical effort. Two general ideas 
have been put forward to try to explain some of the 
observations. One relies on the role of the disorder in 
generating inhomogencity in structurally homogeneous 
samples [15|-[l9j. Within this framework, the crossover 
from insulating to superconducting behavior occurs when 
the Josephson-coupling between superconducting islands 
succeeds in forming a percolation path throughout the 
system. A second scenario invokes the existence of uncor- 
related pre-formed electron pairs which do not constitute 
a condensate but are characterized by an energy gap that 
is associated with the pair binding energy [20[. Other 
models adapt concepts from both pictures and suggest 
that the insulating film is composed of small supercon- 
ducting islands that are uncorrelated and are too small 
to sustain bulk superconductivity [2lJ, |22| . Clearly, addi- 
tional experimental results are required to help shed light 
on this " supercoducting insulator" . 

In this letter we present an experimental study of the 
DOS and corresponding transport characteristics of two 
disordered films on both sides of the transition. These 



measurements show that a similar energy gap exists in 
both the superconducting and the insulating states. We 
present the dependence of this gap on disorder and mag- 
netic field and discuss the possible implications of these 
experimental results towards the understanding of super- 
conductivity in highly disordered films. 

For measuring the DOS of the films wc fabricated tun- 
nel junctions in the following way: A 30 nm Al stripe 
was thermally evaporated on a Si/SiO substrate and was 
allowed to oxidize for a few hours in ambient conditions. 
Subsequently, a 31 nm thick indium oxide (InO) stripe 
was e-gun evaporated perpendicular to the Al stripe, thus 
forming a planar tunnel junction with barrier dimensions 
of 1mm* 1mm. In order to produce InO films with dif- 
ferent disorder, dry oxygen was injected into the evapo- 
ration chamber at different partial oxygen pressure [23 1. 
This resulted in amorphous yet structurally homogeneous 
films with different degrees of disorder characterized by 
their sheet resistance Rn ■ 

The results presented in this letter were obtained on 
two amorphous InO films; one exhibiting superconduct- 
ing transport (sample S) and the other exhibiting insu- 
lating behavior (sample I). Fig. 1 shows the resistance 
versus temperature of both films. It is seen that while 
sample S shows a clear superconducting transition with 
Tc ~ SK, the resistance of sample I increases rapidly 
with lowering T down to our base temperature of 25mK, 
thus showing insulating properties. A fit to Arrhenius 
law: 



Roexp{To/T). 



(1) 



is shown in Fig 1. From this fit one yields an activation 
temperature To = OAK [23|, and Rq = ll.bkfl which 
is larger than the quantum resistance for pairs, 1/Gq = 
h/2e^ ^ 6A5kn. 

Despite this clear difference in the transport proper- 
ties, the tunneling spectra of the two films appears to 
be surprisingly similar. Fig. 2 depicts dl/dV versus V 
curves for both samples. All tunneling measurements 
presented here were performed by standard lock-in tech- 
niques while making sure that the junction resistance was 
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FIG. 1: (color online). Top: Temperature dependence of 
Ra for sample S. The corresponding AlO barrier resistance 
is IMfi. Bottom; Temperature dependence of Ra for sample 
I (heavy red line) and a fit to Eq. 1 (light black line) from 
which To and _Ro are extracted. The corresponding AlO bar- 
rier resistance is 0.7Mf2. Inset: ln{R) versus 1/T for sample 
I. 

at least an order of magnitude larger than the InO sheet 
resistance so that the film could be regarded as an equi- 
potential electrode. Since these are disordered films one 
has to take into account that electronic interactions cause 
a suppression of the density of states (DOS) at low en- 
ergies either due to the Altsuler Aronov [25| zero bias 
anomaly (ZBA) mechanism for weak disorder or due to 
the coulomb gap [l^ for high disordered systems. In 
order to isolate the superconducting contribution to the 
DOS we normalized the curves by the tunneling spectra 
of the films taken at a magnetic field, H, of IIT. The 
justification for this procedure relies on the assumptions 
that superconductivity is fully suppressed at H ~ llT, 
and that the magnetic field hardly affects the normal 
state DOS. Indeed, the dl/dV-V curves of both samples 
exhibit very small magnetic field dependence at energies 
larger than the superconductive gap, and for H > 7.5T 
(where superconductivity is suppressed) the curves are 
practically indistinguishable. 

The dl/dV-V curves shown in Fig. 2 demonstrate that 
a superconducting gap exists in both samples, despite 
the fact that sample I clearly shows insulating transport 
behavior. In order to extract the value of the supercon- 
ducting gap. A, we fit these curves to the BSC expression 
modified by a broadening parameter F that accounts for 
the finite scattering time of the superconducting quasi- 
particles [27[: 



NsiE)/NM{E) =^[{E- iT)/[{E - zVf - A^]'/'-} 

(2) 

The fits depicted in Fig. 2 show that the experimental 
curves of both films deviate from the BSC theoretical pre- 
dictions. In particular, the so called "coherence peaks" 
at the gap edges arc suppressed in our disordered films. 



Nevertheless, the best fits yield A = 0.7m.V for both films 
p8| . The only observed differences between the samples 
are a smooth suppression of the coherence peaks and an 
increase of the sub-gap conductance. 
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FIG. 2: (color online). Normalized tunneling density of state 
obtained at IK for sample I (full red circles) and sample S 
(empty blue circles). The solid lines are fits to the BSC ex- 
pression of Eq. 2 for sample I (dashed red line) and sample S 
(solid blue line). Inset: Raw dl/dV versus V data for sample 
I at H=0 (heavy red line) and H=11T (light black line). 



The fact that the superconductive gap persists into 
the insulating side of the SIT was experimentally im- 
plied in the past. Scanning tunneling microscopy mea- 
surements of the local DOS were performed on disordered 
TiN [i, [H and InO superconductors with different 
degrees of disorder. In both cases it was found that A 
does not decrease with disorder as would be expected 
from the decrease of Tc extracted from the transport 
measurements. In the TiN films A decreased with in- 
creasing disorder, however ^ was found to decrease as 
the films approached the SIT. In the InO samples the 
average A showed no clear dependence on disorder. In 
both cases an extrapolation would predict a finite gap in 
the insulator. Our measurements confirm this trend since 
we observe the same gap magnitude on the two sides of 
the transition. This observation is in accordance with a 
recent theoretical work [g^] which predicts that A is ex- 
pected to remain unchanged through the SIT and may 
even grow in magnitude deep in the insulator. 

As may be expected, the magnetic field, H, has a signif- 
icant effect on A. Fig 3 shows the dependence of dl/dV 
versus V curves on H for samples S and I. It is seen that 
for both samples A decreases monotonically with H and 
is wiped out altogether at H ^ 7.5T. It is interesting 
to relate these results to the magnetoresistance of the 
samples. Fig 3 shows that both samples exhibit a mag- 
netoresistance peak similar to that reported in the past 
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FIG. 3: (color online), (a) Normalized tunneling density 
of state obtained at T=0.5K and at different magnetic fields 
for sample I. inset: the corresponding magnetoresistance. (b) 
Normalized tunneling density of state obtained at O.IK and 
at different magnetic fields for sample S. upper inset: The 
corresponding magnetoresistance for sample S. lower inset: 
A versus H for sample I (empty black squares) and sample S 
(full black circles), obtained from a fit to the BCS expression 
of eq.(2). 

[^-Q- The magnetic field at which the peak occurs, Hp, 
is lower in sample I (1.7T at IK compared to 6.5T for 
sample S). This is consistent with the trend seen in su- 
perconducting samples where Hp was found to decrease 
with increasing disorder (29| . Comparing the tunneling 
and the magnetoresistance results reveals the fact that 
while in sample S the gap is wiped out at fields lower 
than Hp, in sample I A persists up to fields that are 
larger than Hp. It appears, therefore, that the DOS is 
unaffected by the magnetoresistance details, and partic- 
ularly by the resistance peak, and only depends on the 
value of H. This demonstrates, again, that while the re- 
sistance of the two samples may behave differently their 
DOS is very similar. It has been suggested [17[ that 
the MR peak is due to the formation of superconduct- 
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FIG. 4: (color online) . Left: Tunneling density of state versus 
V (linear scale (Left) and logarithmic scale (right)) obtained 
at IK and IIT for sample I (full red circles) and for sample 
S (blue line). 

ing islands with sizes that depends on the magnetic field. 
In this picture. Hp corresponds to the magnetic field 
at which it is preferable for the current to flow through 
normal regions rather than through the superconducting 
islands. Our results are consistent with this picture since 
it implies that the resistance would be sensitive to Hp 
but the DOS would not be affected by it. 

The results presented above strongly suggest that the 
disorder driven superconductor to insulator transition is 
of a geometrical percolation nature. The fact that the 
DOS changes smoothly across the transition and that the 
order parameter is hardly affected by it is consistent with 
a model of superconducting islands of which the effective 
sizes gradually decreases with disorder. The observation 
of coherence peaks in the insulator suggests that super- 
conductive grains that sustain a coherent condensate are 
present in this phase. We envision that the average size of 
these islands decreases with increasing disorder so that a 
larger percentage of islands falls below the Anderson limit 
for superconductivity. This causes the coherence peaks 
to shrink with disorder. Further increasing of the disor- 
der is expected to cause the coherence peaks to vanish 
altogether, however, it is very difficult to perform tun- 
neling junction experiments in this regime of high film 
resistance. 

This "granular" model is further supported by the re- 
sults of the DOS in the normal state. As noted, for 
H > 7.5T the superconductive gap is fully suppressed. 
At these fields dl/dV is found to be proportional to ln(V) 
(see Fig. 4). This behavior is consistent with the zero 
bias anomaly model for weak disorder Due to the 

nearby Al electrode the long range Coulomb interactions 
are screened. Hence for low voltage and temperature, 
such that the thickness of the sample is smaller than 
the thermal length, the DOS is expected to follow the 
screened 2D expression: 
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(00) = -^ In — In- (3) 

where e = max{V,T), z^(oo) is the DOS at high en- 
ergies, T is the inelastic relaxation time, k is the dielec- 
tric function, n is the inverse screening length, b is the 
barrier thickness and g is the dimensional conductance 
G/Gq. Since the slopes of the dl/dV versus ln(V) curves 
are proportional to 1/g, our results enable us to extract 
the conductance ratio between the two films. This yields 
^ = 1.6. We stress that this ratio docs not depend 
on the exact model for ZBA as long as the 1/g depen- 
dence is valid. On the other hand, the conductance ratio 
extracted from the transport at H=11T gives ^ = 11. 
This dramatic difference can be interpreted as an indica- 
tion for granularity in the film. While the transport g is 
sensitive to the current carrying network, tunneling may 
take place into the metallic regions which have a rela- 
tively high conductance even in the insulating phase. In 
this picture the difference between the DOS in the metal 
and in the insulator may be very small. 

In conclusion. We have provided direct evidence for the 
existence of a superconductive gap in an insulator state. 
The gap amplitude is hardly affected by the dramatic 
change in the global transport properties. The results 
strongly suggest that homogeneously disordered super- 
conducting films near the SIT contain superconducting 
islands with bulk properties that are similar on both sides 
of the transition. The DOS is not influenced by the SIT 
nor by the magnetoresistancc peak indicating that the 
insulator contains regions of finite order parameter am- 
plitude that arc similar to those in the superconductor. 
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